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The magneto-optical 
activity of superconducting YBa2Cu307 observed by Karrai 
et al. is not present in many commonly employed models 
of vortex dynamics. Here we propose a simple, unifying pic- 
ture for the frequency dependent magneto-optic response of 
type-II superconductors at low temperatures. We bring to- 
gether Kohn's theorem, vortex core excitations, and vortex 
pinning and damping into a single expression for the conduc- 
tivity tensor. The theory describes magneto-optical activity 
observed in infrared transmission measurements of thin films 
of YBaaCusOr. 

PACS numbers: 74.25.Nf, 74.60. Ge, 74.25.Gz, 76.40. +b 



For applications of superconductivity vortex motion 
and pinning is of great importance. Many models for 
vortex dynamics have been developed such as those of 
Bardeen and Stephen 0], Nozieres and Vinen [||, and 
Clem and Coffey §. One way of studying vortex dy- 
namics is to go to frequencies that are large compared 
with the characteristic frequencies associated with pin- 
ning and damping to observe the free inertial response of 
the vortices. This regime is not well studied and more- 
over a microscopic approach is needed. We note that, 
for high temperature superconductors, estimates of 
the frequency range over which one expects electromag- 
netic absorption (10^° — 10^'' s~^) overlaps with the mi- 
croscopic energy scale associated with quantized quasi- 
particle states in vortex cores [D . Because of the short 
coherence length in high temperature superconductors, 
the energy scale of quasi-particle states localized at the 
vortex core is large. 

These considerations have prompted a re-examination 
of vortex core states. Karrai et al. examined the mag- 
netic field and frequency dependent infrared transmis- 
sion coefhcient of thin films of YBa2Cu307 and found 
evidence for dipolc transitions in vortex cores In 
addition, magneto-optical activity of the superconduct- 
ing state was observed in these experiments . For fre- 
quencies above the vortex resonance this chiral response 
is consistent with the cyclotron resonance of the mixed 
state. According to Kohn's theorem for an isotropic, 
homogeneous electron system in a uniform applied mag- 
netic field H, the only excitation produced by a uniform 
electrodynamic field is the cyclotron resonance at the fre- 
quency ujc = eH/mc, where m is the bare band mass. 



Thus, for an ideal, pure superconductor, cyclotron reso- 
nance is expected. From these considerations, one of the 
requirements of a theory of vortex dynamics is that it be 
consistent with Kohn's theorem. Commonly employed 
theories |0,||j4| fail this test since they are non-chiral. 

The matrix elements and selection rules for the dipole 
transitions in vortex cores have been studied by Janko 
and Shore |^ and Zhu, Zhang and Drew jl^. Kopnin 
proposed a peak in the electromagnetic absorption 
of pure superconductors due to these transitions. How- 
ever Hsu showed that these intravortex transi- 
tions are not excited by long wavelength probes in the 
clean limit. Instead the collective cyclotron motion of 
the center of mass is excited. The {q = 0) conductivity 
is cTxx — {ne^ /m){5{uj — ojc) + 5(ijJ + uJc))/2. The cy- 
clotron resonance exhausts the sum rule dwRe axx — 
(7r/2)ne^/m so that there is no spectral weight at the 
bare dipole transition frequency in the vortex core, 



c{Hc2/ H), which is distinct from luc- There- 



fore, this theory is consistent with the Kohn's theorem. 

According to Hsu |12 1^ pinning activates the oth- 
erwise invisible vortex core excitation. Pinning also 
produces a zero frequency delta function in Re axx- 
These features result from breaking translation invari- 
ance, which invalidates Kohn's theorem. Consistent 
with this picture is the observation that films of 
YBa2Cu307 grown on LaAlOa substrates have a much 
smaller signal at the vortex resonance than those grown 
on Si (with a YSZ buffer). Si/YSZ is a poorer lattice 
match, and more conducive to defect formation. 

Our conductivity is based on identifying coherent ex- 
cited states of a vortex core with translation of the core 
|]l3| . With this identification we relate certain low-energy 
excitations with the velocity v^, of the vortex core. The 
calculation of the q — conductivity is based on three 
equations derived for the clean limit 



J = ?ie(vs + <i>(vL - V5)) 



(1) 



e£ 



vs - Ny{h/2e)wL x z = vg - cj^vl x z , (2) 



vl = vs - (1 - $)Oo(vL - V5) X z 



(3) 



J is the spatially averaged current density, n is the carrier 
density, V5 is the uniform background superfluid veloc- 
ity, = A^/Ep ~ ?i^/2m^^ is the spacing between 
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quasiparticle levels in the core (i5| , A is the bulk gap en- 
ergy, Ep is the Fermi energy, ^ is the coherence length, 
<f> = ujc/^o ~ H/Hc2 is the magnetic field, ujc is the 
cyclotron frequency, £ is the spatially averaged electric 
field, Ny is the areal density of vortices, Ty is a vortex 
damping rate, and a is a pinning frequency. 

Eq. (|^) is obtained by evaluating the expectation value 
of the current operator for a vortex core of velocity 
and then applying Galilean invariance. This ensures that 
when Vi = V5 the current is simply that of all carri- 
ers moving at V5. The — V5 part of the current, 
not present in conventional theories, comes from coher- 
ent excitations of the core jl^ . In our calculation of the 
current due to vortex core excitations a backflow is not 
included when vl 7^ V5 since V • J(r) = —dp/dt ^ 
for the a.c. case. This differs from the usual low fre- 
quency, quasi-static, vortex motion calculation. In our 
calculation we have focused on the polarizability of the 
core which couples to an applied electric field. We have 
neglected modifications due to screening. 

Eq. (H) says that the total electric field is the sum of a 
gauge term from the uniform acceleration of the super- 
fluid background plus the average Josephson electric field 
due to transverse motion of the vortices. The sign of the 
latter term corresponds to a magnetic field pointing in 
the +z direction and positive charge carriers. 

Eq. (|^) is an equation of motion for a vortex core. It 
was derived in Rcf. p^ ] by looking at the microscopic 
equation of motion for low energy quasiparticles in the 
core and applying the gap equation. Ty is identified as the 
quasiparticle relaxation time in the vortex core. It is the 
analogue of the Drude relaxation time for itinerant elec- 
trons. Note that the damping in Eq. (^) is inversely pro- 
portional to Ty whereas in theories based on the Bardeen- 
Stephen model 0^ it is proportional to the normal state 
electron transport lifetime r. This difference is due to the 
clean limit {QoTy ^ 1) in our theory and the effectively 
dirty limit {Qqt <^ 1) conditions of the hydrodynamic 
models. The Bardeen-Stephen model breaks down when 
r ^ 00. Nevertheless, as we note below, the two ap- 
proaches agree in certain cases. The harmonic pinning 
term, — a^r, can be derived from a simple single-particle 
short-range repulsive potential [p^ . It has an appealing 
and generally expected form (for low vortex densities). 
Notice that in the limit of small dissipation and pinning 
the steady state solution of Eq. (||) is v^, = V5 as required 
by Galilean invariance. 

Eq. (|^) contains a term corresponding to a Magnus 
force. The existence of a Magnus force has been contro- 
versial in type II superconductors Recently, however, 
very general arguments have been given 1 17| for a Magnus 
force even in the presence of pinning and viscous drag. A 
comparison jl^ of Eq. (||) to the Magnus force gives an ef- 
fective inertial "mass" of the core Aly = hn/2ilo « n^'^m 
per unit length in the low vortex density limit. This 
mass, determined microscopically from the core energy 



level spacing, is not equal to the Suhl value ||18|,^ for the 
vortex mass. Vortex dynamics studies by Kopnin in su- 
perconductors 1 11] and superfluid '^He and by Baym 
and Chandler [21| in superfluid ^He find a vortex mass 
similar to ours. The different mass obtained by Suhl ap- 
pears to be a consequence of the neglect of quantized core 
excitations in the Landau-Ginzburg theory. The electro- 
dynamic response we are considering corresponds to a 
polarization of the vortex core and not the motion of a 
current pattern across an otherwise uniform charge den- 
sity as discussed by Suhl. 

The factor (1 — $) in Eq. (j^) was not derived micro- 
scopically but included in order that Re (Txx{tjj) > [p^ . 
This is ensured because Eqs. (^) and (||) lead to 



2"'dt 



(1-*)|V5 P] 
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nm<E> 



(4) 



Taking the time average over one period shows that 
(3 ■ £) > 0. We note that other plausible generaliza- 
tions of conventional vortex dynamics theories to include 
the vortex mass generally do not satisfy this condition. A 
decreasing ilp with magnetic field is expected on general 
grounds. It can come from a renormalization of the Mag- 
nus force due to the presence of nearby vortices. Also the 
magnetic field increases as the vortices overlap and will 
affect the dynamics. These effects were ignored in our 
single vortex calculation. 

The three equations in four unknowns (J, £, vl, V5) 
may be used to eliminate vl and V5 leaving J = a£ with 
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uj{uj T ^0) + (1 - ^){iuj/Ty - a^) 

(UJ =F flo){uj ± LOc) + iuj/Ty — oP- 



(5) 



where g± = a^x i iu^y and Hq = (1 — The ± 

refers to cyclotron resonance active (— ) or inactive 
modes of circularly polarized light. This conductivity 
satisfies a number of important limits. First, when B — Q 
the London conductivity is recovered. For zero pinning 
and zero dissipation the conductivity contains simply the 
cyclotron resonance as given above. In the limit of no 
pinning, zero frequency, and low vortex density, pxx = 
(m/ne^ri,)<I> = p^J^ (p"^. is the normal state resistance) 
since Ty is identified with r. Therefore we obtain the 
same expression in the clean limit as was found in the 
hydrodynamic theories In this limit we also obtain 

Pxy — B/nec as expected from Galilean invariance. The 
sum rule is satisfied. For example pinning produces a 
delta function, 



ne2 [ (1 - $)a2 

Re cTxx[uj) = 

m 



TrS{uj) . 



(6) 



This spectral weight is removed from the cyclotron reso- 
nance. The cyclotron resonance is shifted to lUc + ct'^/^o 
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to produce a hybrid cyclotron-pinning resonance. Spec- 
tral weight is also transfered between the cyclotron res- 
onance and the vortex core resonance. In Fig. || we plot 
the conductivity and show how it is affected by pinning. 
For values of the pinning frequency that are consistent 
with experiments our theory predicts a very small opti- 
cal response at in the cyclotron inactive mode. 

If we calculate the thin film surface impedance, in 
the low frequency limit (lo <C a, ^o) and the limit where 
pinning dominates the magnus force (o? ^ t^f^o) we ob- 
tain the Gittleman-Rosenblum , Coffey-Clem 1 23 1 re- 
sult for pinned vortices. Alternatively we could derive 
this by dropping the vl, v^^ and x z terms in Eq. 
(^. Dropping the vl term is equivalent to setting the 
mass to zero. Using Z ~ l/o'±(aj)t, [t is the film thick- 
ness), and defining a drag coefficient r\ = M^/ry, and a 
pinning force constant k = M^a^, 



Z 



(7) 



The first term is the London reactance term. It differs 
from the Coffey-Clem result by the 1 — $ factor. This re- 
duction of the effective superfluid density by 1 — {H/ Hc2) 
is not unexpected nor contradicted by experiment. Eq. 

was obtained under the clean limit assumptions of 
our theory. Remarkably it has the same form as the dirty 
limit case. Therefore if the Magnus force were present in 
the dirty limit, as argued by Ao and Thoulcss then 
the correct conductivity function in that limit might also 
resemble Eq. (||). 

We have fit measurements of the infrared transmis- 
sion of thin films of YBa2Cu307 with this theoretical 
conductivity. The experimental technique is described 
in Refs. P,|7|,p^ . The transmission coefficient is = 
47V/ 1 TV + 1 + {Att / c)a±-0? . N is the refractive index of 
the substrate. In Fig. | we plot the ratio of the trans- 
mission coefficients T"*" and T~ . This ratio mostly can- 
cels out non-chiral components of the transmission. The 
figure shows a definite optical activity in the sample 
which at high frequencies is well described by a lossless 
free electron conductivity [0,0. At frequencies below 
^ 50 cm~^, however, the chiral response cannot be well 
represented even if damping is included in the free elec- 
tron conductivity as can be seen in Fig. In particular 
there is a steep drop in /T^ at low frequencies which 
corresponds to the onset of the hybrid cyclotron-pinning 
resonance. This feature is well fit by our theoretical con- 
ductivity. The free parameters in our fits are t^, a, and 
f2o- We note that the conventional conductivity func- 
tion [^^, even when generalized to finite vortex mass, is 
non-chiral and would predict simply /T^ = 1. That 
/T^ falls below unity at 30 cm~^ is particularly note- 
worthy as it implies a sign reversal of the a.c. Hall effect 
at low frequencies. 

Another observation from the data displayed in Fig. |^ 
is that the chiral signal associated with the vortex core 



resonance is very small compared with the cyclotron- 
pinning resonance. Our theory also predicts a very small 
chiral signal but does not agree in detail with the mea- 
sured signal. The measured resonance is larger than the 
theory predicts. Moreover, there are substantial non- 
chiral features observed in the transmission spectrum 
around 65 cm~^ . These features are brought out in unpo- 
larized transmission as can be seen in Fig. H The peak 
around 65 cm^^ is substantially larger in the unpolar- 
ized spectrum than in the T'^ /T~ spectrum. This figure 
also shows the theoretical unpolarized transmission. The 
corresponding small broad peak in the calculated curve 
is chiral and has been adjusted to have the maximum 
amplitude consistent with the low frequency rise in the 
transmission. Our conductivity does not contain the non- 
chiral response observed in the experiments. Neverthe- 
less, we believe jlj] that this non-chiral response and the 
enhanced chiral response may be due to vortex core exci- 
tations. The effects of defects beyond harmonic pinning 
induce optical transitions which, because of the break- 
down of the cylindrical symmetry, have less restrictive 
optical selection rules. 

The presence or absence of polarization dependence 
is controversial. The considerations of Zhu, Zhang and 
Drew and Janko and Shore would predict that 
the vortex resonance is in the cyclotron active mode. The 
conductivity here leads to the opposite effect: a weak vor- 
tex resonance in the cyclotron resonance inactive mode. 

In conclusion, we have presented a conductivity func- 
tion which produces magneto-optical activity at high fre- 
quencies and explains how it is diminished by vortex 
pinning and damping at low frequencies. The theory 
predicts that the dipole-allowed vortex resonance is very 
weak due to strong screening by the vortex motion even 
in the presence of pinning. Thus, the Kohn theorem is 
very robust in this system. The conventional theories 
of vortex dynamics are not consistent with Kohn's the- 
orem and do not contain optical activity. Therefore our 
conductivity function provides a better starting point for 
studying the interaction of vortices with impurities and 
the excitation of vortex core states. 
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FIG. 2. The chiral response for the YBa2Cu307/Si sample 
at J? = ±12 T and 2.2 K as determined from measurements 
with a polarizer with a 0.9 mm quartz wave plate. The solid 
curve is the best fit to our conductivity function. The best fit 
parameters are a = 50±2 cm^^, l/r„ = 40±5 cm~^, when fio 
was taken as 60 cm~^ and the cyclotron mass m was taken to 
be 3.1 electron masses [7]. The dashed curve is the best fit to 
the damped free electron model. The corresponding damping 
paramter is 1/t = 56 cm"'^. 

FIG. 3. Transmission ratio at 2.2 K for the YBa2Cu307/Si 
sample taken at 12 T and 2.2 K in unpolarized light. The solid 
line is the theoretical result using the best fit parameters from 
the analysis of the chiral data in Fig. 2. 



FIG. 1. The frequency dependence of the conductivity 
function, Eq. (5), is plotted for cyclotron inactive (-I-) and 
active ( — ) modes. The parameters arc, loc ~ 10 cm^^, 
a = 20 cm~^, Q.o = 60 cm~^, and 1/t„ = 10 cm~^. There 
is also a delta function at zero frequency whose strength is 
given by Eq. (6). 
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